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Abstract 

This paper studies the effects of axial velocity profiles on vortex decay, in order to properly initialize and 
simulate three-dimensional wake vortex flow. Analytical relationships are obtained based on a single vortex 
model and computational simulations are performed for a rather practical vortex wake, which show that the 
single vortex analytical relations can still be applicable at certain streamwise sections of three-dimensional 
wake vortices. 


Introduction 

Aircraft trailing vortex wakes can cause seri- 
ous loss of control when following aircraft encounter 
them. The following aircraft can be subjected to 
rolling moments which exceed the aircraft’s roll con- 
trol authority, leading to a dangerous loss of altitude, 
and to possible structural failure. With the advent 
of large transport aircraft, the wake vortex problem 
has taken on added significance. Since the probabil- 
ity of an aircraft-vortex encounter is greatest in air- 
ports where aircraft operate in close proximity, air- 
craft spacing at congested airports is dictated by the 
characteristics of the vortex wakes left in the terminal 
area 6 . Therefore, prediction of wake vortex trajecto- 
ries and strengths is especially important for effective 
airport flight control and maintain maximum traffic 
volume. Both government agencies (NASA, FAA) 
and industry (Boeing) feel an urgent need to solve 
the problem. 

In order to fully assess the hazard, complete in- 
formation is needed regarding such processes as the 
initial organization of the wake vortex sheet, the sub- 
sequent roll-up and generation of primary vortices, 
the descent and decay of the wake in a real atmo- 
sphere, and the response/control characteristics of 
the encountering aircraft. Since wake vortex trajec- 
tories and strengths are altered radically by interac- 
tions with the ground plane and by atmospheric con- 
ditions, computational simulations have been focused 
on the viscous interaction between vortex wakes and 
the ground plane, including atmospheric effects such 
as stratification, wind shear and turbulence. An in- 
depth understanding of the mechanisms that bring 
about wake transport and decay can provide infor- 
mation which could enable the alleviation of airport 


congestion. 

Previously, most computational simulations have 
been two-dimensional or quasi-two-dimensional (e. 
g., Ref. 2, 4, 10, 11, 13-17). Those simulations are 
applicable when the streamwise vortex axial veloc- 
ity gradients are small comparing with the velocity 
gradients in vertical and spanwise directions. Two- 
dimensional simulations can predict the vortex trans- 
port behavior at a satisfied level in comparison with 
measurement data. However, wake vortex decay his- 
tories calculated in two-dimensional simulations are 
questionable. Phenomena such as vortex stretching, 
vortex bursting or breakdown and Crow instability, 
which are intrinsically three-dimensional, are unable 
to be captured using two-dimensional models. 

The study will then be concentrated on three- 
dimensional simulations. Since it is an unsteady 
problem, flow initialization is required. To spec- 
ify a physically sound initial three-dimensional wake 
vortex flow field is always a challenging problem in 
fluid dynamics, since the roll-up process itself is a 
complicated phenomenon. However, since only far- 
downstream flow field is of interest in wake vortex 
hazard predictions, detailed information of how the 
vortex is generated on the wing is another research 
topic and thus excluded in this research. Rather, the 
initial vortex flow field will be analyzed using integral 
relationships, instead of detailed complicated CFD 
approach, although and latter and some experiment 
data are probably necessary to evaluate the integral 
model and calibrate some coefficients. 

In the three-dimensional algorithms, the up- 
stream in-flow boundary conditions are specified as 
the initial vortex flow field and thus again it shows the 


‘Assistant Professor, Mechanical Engineering Department, AIAA Member, 

Copyright ©1997 by the American Institate of Aeronautics and Astronautics, Inc. All rights reserved. 


1 




importance of a well-defined initialization. We ran a 
simulation with uniform axial velocity component (no 
axial shear) to test the development of the axial flows. 
Figure 1 is a two-dimensional section of axial veloc- 
ity contours, with x the axial direction, y and z the 
lateral and vertical directions, respectively. It shows 
that a non-physical axial flow field, with half of the 
axial flow towards the airplane (negative value) and 
half of it away from the airplane (positive value), has 
been generated. That means non-practical axial flow 
initialization can result in incorrect three-dimensional 
simulations. 

Initially, a completely rolled-up wake vortex pah- 
can be modeled as two non-interacting counter- 
rotating vortices, because of the large separation dis- 
tance between the two vortices (about eighty percent 
of aircraft wing span) and the small vortex core ra- 
dius (about five percent of the wing span). Several 
azimuthal velocity models have been studied (such 
as Ref. 7) and most of them are agreeable with mea- 
surement data if the model coefficients are properly 
adjusted. However, the axial velocity component has 
been of some difficulties to determine, due to lack of 
systematical measurement data and reliable compu- 
tational simulations. 


profiles are represented in Figure 2. All three types 
have been observed in flight tests. 

Theoretical Analysis 

The importance of the axial velocity component 
in three-dimensional simulation can be explained by 
a perturbation expression based on two-dimensional 
formulations, which can be considered as a three- 
dimensional case with uniform axial velocity. In two- 
dimensional simulations, the axial direction is con- 
verted to time with t = xjU 0 , where U a is the 
flight speed of the aircraft and x is conventionally 
assigned to axial direction and y and z are span- 
wise and vertical directions, respectively. Hence, a 
two-dimensional calculation is actually a steady state 
three-dimensional calculation with a uniform axial 
velocity component, which is the flight speed in this 
case. 

Consider the transport equation for the axial com- 
ponent of vorticity, ( x , which has the direct contribu- 
tion to the vortex roll hazard. Strong axial currents 
occur near the vortex core, where Cr concentrates. If 
we put a small perturbation to the uniform axial ve- 
locity, the axial vorticity transport equation becomes: 


There are several researchers who have studied 
vortex axial velocity after completion of vortex roll- 
up. Batchelor 3 analyzed the dynamical necessity of 
axial flow for inviscid vortices: the radial pressure 
gradient balances the centrifugal force, where any 
change in the azimuthal motion with distance down- 
stream produces an axial pressure gradient and con- 
sequently an axial acceleration. Moore 8 brought vis- 
cous effects into the analysis and estimated the axial 
flow deficit due to the boundary layer on the wing. 
He concluded that both jet and wake type axial flows 
could occur depending on the distribution of loading 
on the wing. Brown 5 showed that axial flow effects of 
the wing profile drag and lifting system depended on 
the ratio of profile drag to induced drag. Recently, 
Rule and Bliss 12 proposed a model based on the Betz 
method. They obtained a set of ODEs to calculate 
both tangential and axial velocity profiles. 
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where 

Cr =dW/dy-dVfdz, 

C y = du/dz - dw/dx 

and 

Cr = dVjdx - du/dy . 

Substituting the expressions of Cy and C* suid making 
use of the continuity 


du dV dW 
dx + dy ^ dz 


Eq. (1) becomes: 


( 2 ) 


In spite of the complex of initial axial flows, at the 
near field (close to the flying aircraft), the profiles of 
the axial flows can be categorized into three types 
(in coordinates attached to the aircraft): (1) jet type 
(excessive to the free stream); (2) wake type (deficit 
to the free stream); ^3) combinations of the two (with 
jet type close to the vortex center and wake type fur- 
ther away from the center in radial direction). Those 
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Divide the above equation with the constant flight 
speed and define u' = u/U 0 , we have 
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1 >dx + u c ( ay + 8z > 

dW du' dVdu' v „ 0 

dx dy + + • (4) 

Equation (4) is an exact equation for axial vorticity 
component. For analytical study purpose-, here we 
assume that the axial flow deviation is small compar- 
ing with the flight speed, i.e., u' < 1. The equation 
that we obtain is 


vortex velocity distribution. Since most of the axial 
velocity deviation is inside the vortex core, both u' 
and 8u' /dr are approximated to be zero outside the 
vortex core. Substituting Eq. (8) into Eq. (7) we 
have 
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After integration by parts and re-arrangement, we 
obtain 
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It can be seen that the relationship between total 
circulation decay and the vortex core size change is 
v. / depending on the axial velocity inside the core. 


Now consider a single vortex case. We integrate Eq. 
(5) throughout each y — z section, noticing that the 
area integration of axial vorticity is the total circu- 
lation at that section. It can be shown that the re- 
sultant equation for the total circulation (at certain 
x location) is: 


If the axial flow is wake type ( u ' < 0), when the 
circulation decreases, the core radius decreases. If the 
axial flow is jet type (t/ > 0), the numerator in the 
fraction at the right-hand side of Eq. (10) is nega- 
tive, while the denominator is still positive for u' < 1. 
Thus, when the circulation decreases, the core radius 
increases. 


<*Too_ f , dWdu' , dVdu\ 1M 
dx 7/ dx dy + dxlh' dA ' 


( 6 ) 


All the integrations of the convection terms (in con- 
servative forms) and diffusion terms are zero, due to 
the fact that the axial vorticity and its derivatives 
are zero at the infinite boundary. We also note that 
Eq. (6) reduces to Betz’s first invariant 7 if the axial 
velocity component is uniform in y — z plane and x is 
replaced by tfU 0 . That is, for two-dimensional single 
vortex cases, dT^/dt — 0. 


For the single vortex case, axi-symmetry can be 
assumed. Using cylindrical coordinates, Eq. (6) be- 
comes: 

_ f 3Ve (r, x) du'(r,x) ,, 

dx J A dx dr dA ’ ^ 

where V 6 is the tangential velocity of the vortex. 
Since we are trying to deduce the qualitative influ- 
ence of the axial velocity profile on the vortex decay 
behavior, here we choose Ranldne vortex to derive a 
theoretical relationship, without losing generality. In 
Rankine vortex, 


Then we look at the relationship between maxi- 
mum tangential velocity and the vortex core under 
the influence of axial velocity. Substitute the expres- 
sion 


Foo 27rr c V0 max (11) 

into Eq. (10), we can have: 
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( 12 ) 


Since u f < 1, the term in the parenthesis at the right- 
hand side of Eq. (12) is always negative. Therefore, 
change of maximum tangential velocity is always in 
the opposite sense of change of vortex core radius, 
regardless what type of the axial velocity profiles (jet 
or wake). When the vortex core grows, the maximum 
tangential velocity decreases, which is the same case 
as in two-dimensional single vortex. 


In Ranldne vortex, the vorticity of the vortex ran 
be written as 



when r < r c , where r c is the radius of the vori 
core. Outside the vortex core, V 6 follows potent 



for r < r c and Cc = 0 for r > r c . 
(13) into Eq. (10), we have 


(13) 

Substituting Eq. 
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Again, note that the denominator in the parenthesis 
at the right-hand side of Eq. (14) is positive for any 
[n'l < 1. That means when the vortex core grows, the 
vorticity inside the vortex core decreases, so does the 
maximum tangential velocity as shown in Eq. (12). 


It should be noted that the above relationship (i.e. 
Eqs. (10), (12) and (14)) is valid under the condition 
that u' < 1 for a single vortex. However, it offers a 
guide line to estimate the influence of axial velocity 
profile. 


Computational Model 

In realistic wake vortex simulation for airport traf- 
fic control purpose, several effects have to be in- 
cluded, such as interactions between a pair of vortices 
in the fully rolled up aircraft wake and the ground 
and the terminal area atmospheric conditions. The 
computational model used in this research is based on 
the NASA TASS (Terminal Area Simulation System) 
computational program. 


u> = ^'maxtxpi-cr 2 /r\) , ( 16 ) 

where + and — is corresponding to jet and wake type 
axial flows and c is a axial flow spreading coefficient 
in radial direction. In the test cases here we chose 

u 'max =0-1 and c = 10. 

As mentioned previously, the computation is for 
unsteady, three-dimensional simulations. If the case 
reaches steady state, then the time march is consid- 
ered as pseudo-time march. The steady-state con- 
vergence criterion used is 10 -3 drop of residue for 
all velocity and pressure (although in most cases the 
steady state is hard to reach). The initial condition 
is started with two single vortices (Eq. (15)) with 
axial velocity (Eq. (16)). That initial condition is 
stacked for about 20 percent of the total length of the 
domain in axial direction. Then the inflow bound- 
ary condition is kept as the same as what initially 
specified (at the upstream end of the domain). Both 
the top boundary and downstream-end boundary are 
open boundaries. The bottom boundary is specified 
as ground boundary and no-slip boundary conditions 
are used. 


The code uses Arakawa 1 C-type grid and is con- 
servative quadratically. Meteorological phenomena, 
such as atmospheric turbulence and stratification ef- 
fects are well embedded in the computational model. 
The Smagorinsky sub-grid scale model is used for 
the large-eddy turbulence simulation. The readers 
are referred to the two references by Proctor 9 - 10 for 
the details of the numerical model. Some features 
are briefly stated here. The equations are in primi- 
tive variables for time dependent, compressible flows. 
The explicit scheme uses 4th-order central differ- 
ence on staggered, uniform grid in space, with time- 
split march on acoustically active terms and non- 
acoustical terms. The time march step is dynami- 
cally calculated for each time march to satisfy stabil- 
ity conditions. 

Initially we put the vortex pair high above the 
ground and the separation between the vortices is 
large. Therefore during early time period, the vortex 
behaves similar to a single vortex. In stead of Rank- 
ine type vortex model, the following vortex model 
is chosen, to be compatible with others in the wake 
vortex simulations 10 : 


The axial velocity profiles are set as: 


The initial circulation and flight speed are 
550m 2 / $ and 61.75 m/s, respectively, which are close 
to a DC-10 vortex wake. The initial r c is chosen as 
7m. The initial height from the ground and sepa- 
ration between the two vortices are 70m and 50m, 
respectively. The grid sizes are 1.75m my — z plane 
and 4m in axial direction. The Smagorinsky turbu- 
lence model is turned off in the current runs, by using 
a very small coefficient with the term so that the vis- 
cous effect is mostly from the molecular viscosity. 

Two test cases are shown in the following, one 
with a wake type axial velocity profile and the other 
with a jet type. Since the cross wind effects are not 
included, the flow is symmetric in spanwise direction, 
thus only half of the flow field needs to be calculated. 
The domain is 640m in axial (r) direction, 52.5m in 
spanwise (y) direction and 104m in vertical (z) direc- 
tion. The time march is stopped at 0.3 minute and 
steady state solution is not reached, according to the 
convergence criterion specified. That is probably be- 
cause the separation near the ground is a unsteady 
phenomenon, which we do not include in the discus- 
sion in this paper. 

Figure 3 is the x — z plane wake-type axial ve- 
locity contours at y — 25m (about the position of 
the center of the vortex), which shows the propa- 
gation history for the axial velocity profile at 0.1, 
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0.2 and 0.3 minutes. It can be seen that at 0.1 
minute, the profile reaches about 61.75m/* x 6s = 
370.5m. At 0.2 minute, the profile is supposed to 
reach 61.75m/* x 12s = 741 m, which is larger than 
the domain size. Then at 0.3 minute, the flow field 
is more established than at 0.2 minute and we chose 
that time level to analyze the flow field, although the 
steady state solution is not reached as explained. The 
same propagation type is found in the case specified 
with jet-type axial velocity profile. 


The total circulation should be used to verify Eq. 
(10), which is probably more interesting, because the 
sign of the axial velocity deviation changes the sign 
of the circulation/core-size relation. However, since 
the computation includes two vortices as well as the 
ground plane, total circulation for just one vortex is 
difficult to calculate and mostly not of interest for 
practical purpose. 


Conclusion 


In both cases, the maximum vorticity at the vor- 
tex center decreases and the vortex core size in- 
creases, in agreement with Eq. (14). The vortex 
core size is determined from the span- wise velocity (v) 
contours. Since the vortex has a downward motion 
due to interactions between the two vortices in the 
vortex pair, vertical velocity component is not a good 
measure for vortex core size. The distance between 
the zero span- wise velocity contour and the minimum 
contour (the largest negative value) is used as the vor- 
tex core radius. We consider that way of quantifying 
the vortex cores reducing the influences from both the 
mirror vortex and the ground plane. Figure 4 shows 
. e y z section contours of ^-direction vorticity and 
the u-velocity component at the first * plane, which 
is the same for both wake and jet cases. Figure 5 is 
for the wake case at the last s plane and Figure 6 
is for the jet case at the same plane. Note that the 
u contours and x - vorticity contours are not exactly 
at the same a: location, because of the staggered grid 
used m the computation. It can be seen that the 
maximum vorticity level decreases about 15.2% and 
the core size increases about 13.3% in the wake case, 
while in the jet case the maximum vorticity level de- 
creasesabout 17.9% and the core size increases about 
• Vo. In general, those two cases do not show signif- 
icant difference. We expect that if the domain size 

is longer in the axial direction, more difference will 
show. 


The axial velocity component is important for 
three-dimensional wake vortex simulations. Initial 
axial velocity profiles influence the vortex decay be- 
havior. Analytical relationships for small axial ve- 
ocity deviation have been developed in single vortex 
cases. Those results show that the sign of the axial 
velocity deviation can change the relation between 
the total circulation and the vortex core size, but not 
the relations between vorticity, maximum tangential 
velocity, and the core size change. 

. Computational simulations were tested for both 
jet and wake type axial flow, with minimal influence 
from the ground and the mirror vortex in the wake 
vortex pair. Not surprisingly, the computational re- 
sults are in agreement with the trend predicted in 
the analytical studies for the relationship between 
vorticity and core size, although at some locations 
opposite trend is detected with negligible influence 
on the. overall vortex decay. The relationship for the 
total circulation needs to be validated using different 
simulation setup. For future research, the influence 
of strong axial currents, instead of small deviations 
should be studied. ’ 
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Figure 1: 2D Section of Axial Velocity Contours 
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Figure 2: Axial velocity profiles 
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Figure 3: Axial velocity development at: (a) t=0.1 min., (b) t=0.2 min., (c) t=0.3 min. 
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Figure 4: (a) Axial vorticity (b) Span-wise velocity contours at the first axial plane 
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Figure 5: (a) Axial vorticity (b) Span-wise velocity contours at the last axial plane in the wake case 
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Figure 6: (a) Axial vorticity (b) Span-wise velocity contours at the first axial plane in the jet case 
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